Dispersion interactions and reactive collisions of ultracold polar 

molecules 

Svetlana Kotochigova* 
Department of Physics, Temple University, Philadelphia, PA 19122-6082, USA 

Progress in ultracold experiments with polar molecules requires a clear under- 
standing of their interactions and reactivity at ultra-low collisional energies. Two 
important theoretical steps in this process are the characterization of interaction 
potentials between molecules and the modeling of reactive scattering mechanism. 
Here, we report on the ah initio calculation of isotropic and anisotropic van der 
Waals interaction potentials for polar KRb and RbCs colliding with each other or 
with ultracold atoms. Based on these potentials and two short-range scattering pa- 
rameters we then develop a single-channel scattering model with flexible boundary 
conditions. Our calculations show that at low temperatures (and in absence of an 
external electric field) the reaction rates between molecules or molecules with atoms 
have a resonant character as a function of the short-range parameters. We also find 
that both the isotropic and anisotropic van der Waals coefficients have significant 
contributions from dipole coupling to excited electronic states. Their values can dif- 
fer dramatically from those solely obtained from the permanent dipole moment. A 
comparison with recently obtained reaction rates of fermionic ^"K^'^Rb shows that 
the experimental data can not be explained by a model where the short-range scatter- 
ing parameters are independent of the relative orbital angular momentum or partial 
wave. 

I. INTRODUCTION 

Recent advances in creating ultracold polar molecules in the lowest rovibrational ground 
state [H 12] have created a new scientific playground for studying quantum phenomena that 
govern collisions and interactions between molecules at low temperatures. Some key theo- 
retical insights into the dipolar character of these interactions have already been developed 
over the last few years [3H9]. It was shown that dipole-dipole forces can be quite strong 
and may give rise to complex many-body physics in cold polar gases. Molecules with dipole 
moments, which interact via very long-range dipolar forces when oriented in an external 
electric field, can form strongly correlated condensed matter-like systems, realize effective 
spin models, and create field-linked states [9HT2]. 

Currently, there is also significant interest in investigations of scattering properties of 
ultracold molecules. It is expected that such properties will be very different when molecules 
are held in optical dipole traps or when held in individual sites of an optical lattice. Holding 
molecules in individual sites of an optical lattice prevents them from approaching each other 
and allow dipole-dipole forces to play a dominant role in interactions between polar molecules 
[T3l m]. Intriguingly, in lower dimensional systems, where tight confinement along one or 
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more spatial directions is applied, dipole forces can be engineered to be non-destructive or 
repulsive and collective effects are predicted p^HTS] . For example, Ref. [T7j has suggested 
a technique that decreases inelastic coUisional rates and enhances elastic collisional rates 
using a combination of external electric and microwave fields in a two-dimensional lattice. 

In a three-dimensional trapping environment molecules can scatter freely. Generally, 
polar molecules are expected to be very fragile to the destructive nature of these collisions. 
There is experimental evidence of the significant role that such collisions play in defining 
molecular lifetimes jl9l423] . The physical origin of this loss are rovibrational relaxations 
and reactive collisions, both of which occur at short-range, when the molecules are close 
together. An understanding and quantitative description of these collisions might help 
define conditions under which an ultracold molecular system is long lived. Furthermore, 
learning about the short-range region will help to unravel the role of reactivity in ultracold 
collisions. 

The effect of reactivity on the molecular lifetime at ultracold temperatures, with a few 
exceptions [231 El] , remains largely unexplored. References [251 - I28] are devoted to theoret- 
ical developments towards a full quantum dynamics calculation of reactive collisions with 
cold molecules. The quantum mechanical description of reactions is challenging due to the 
complex nature of the differential equations and boundary conditions as well as difficulties 
with generating potential surfaces of sufficient quality for three and four-atomic systems. 
The most successful theories are associated with light atom-diatom systems [2SHSI] that are 
of astrophysical interest, can be Stark decelerated, or cooled with a buffer gas. Recently, ex- 
tensive quantum scattering calculations in Refs. [321 133] have shown significant influence 
of so-called 'virtual states' in the entrance channel of the collision complex. 

The interaction of molecules without an external electric field is dominated by the dis- 
persive van der Waals forces rather than dipole-dipole interactions. At low temperatures 
quantum mechanical effects play a prominent role in the molecular scattering from such 
potentials and are crucial for the description of the interplay between inelastic and elas- 
tic collisional rates. The present study explores these dispersion interactions, which act at 
short- to medium-range, and introduces a scattering model of ultracold polar alkali-metal 
molecules. We use a characteristic length scale to distinguish between the dispersion and 
short-range interaction regions. The short-range boundary, which defines the lower limit of 
the scattering model, is described by Cq"/R^j. ^ 2Be, where Cq°/R^j. is the isotropic disper- 
sion potential and Be is the rotational constant of the ground state molecule. The short-range 
separation i?sr must be smaller than the isotropic van der Waals length = {2^Cq° /h'^Y^'^, 
where fj, is the reduced mass of the molecules. An accurate estimate of these characteristic 
length scales for KRb and RbCs will be given below. 

Our scattering model is an extension of the theory developed by Mies and Julienne pHl35]. 
Here we modified this scattering theory by making the short-range boundary conditions more 
flexible by assuming that not all molecules that penetrate to short-range will be lost. This 
allows us to introduce two short-range parameters linked to the rovibrational structure of exit 
reaction channels. Recently, a model based on quantum defect theory has been introduced 
in Ref. [2l|, which accounts for the short-range interactions by connecting a complex valued 
or optical potential to the van der Waals potential. 

Since our scattering theory is solely based on the dispersion interaction potentials, we first 
calculate molecule-molecule and atom-molecule van der Waals coefficients by integrating 
products of the dynamic polarizabily over imaginary frequencies. Unlike for atom-atom 
interactions molecular interactions depend on their rotational and vibrational state and can 
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have contributions from transitions within the ground state potential as well as from the 
electronically excited spectrum. For polar molecules the former contribution is nonzero, 
even for levels with a well defined angular momentum while the excited-state contribution 
can be significant. The excited contribution is often missing in the analyses of short-range 
molecular interactions [T71 ES] • 

The paper is organized as follows. In section |TT] we describe the theory of dispersion 
interaction isotropic and anisotropic potentials between molecules and molecules with atoms. 
The numerical potentials and van der Waals coefficients are obtained for the X ground 
state of KRb and RbCs as a function of its vibrational quantum number and presented in 



Section |ITT} Section |IV] gives basics of our scattering model and provide examples of the 
inelastic scattering rates for ultracold colliding molecules. The last Section |V] summarizes 
essential conclusions of the presented research. 



II. DISPERSION POTENTIALS FOR THE MOLECULAR INTERACTION 

We describe the dispersion interaction potential between molecules A and B, each in 
a rovibrational level |X, w JM) = \i,M) of their X electronic ground state, by assuming 
that the molecules are far apart and their wavefunctions do not overlap. Here, the magnetic 
quantum number M is the projection along a laboratory fixed coordinate system of their total 
angular momentum J and i describes all other quantum labels. Their energy Ei only depends 
on the label i and not on M. These assumptions allow us to use a (degenerate) second- 
order perturbation theory similar to that described by Ref. p7j. For two molecules the 
matrix elements of the dispersion potential between the product states \iA,MA',iB,MB) = 
\iA, MA)\iB, Mb) and \iA, Mj^^is, M'^) with different projection quantum numbers but the 
same angular momenta J4 and Jb are 

TT V- {iA,MA;iB,MB\Vdd\SA]SB){SA;SB\Vdd\iAM'^;iB,M'j^) 

tVdispl^j = - 2^ r= — — — — , (Ij 

Sat^^aMa + ^Sb - y^iA + ^iB) 

SbJ^^^B'^^'b 

where the sums Sa and 5*^ are over all electronic, rovibrational, continuum states of 
molecules A and B, respectively, restricted to states with energy Es^ + Esg 7^ Ei^ + Eig. 
The operator Vdd is the dipole-dipole interaction Hamiltonian 



V^(R) = -v^ E ( I i I ) diA^ . (2) 

where d'^ and are the two rank-1 spherical dipole operators of the molecules and R is 
the separation between and orientation of the two molecules with respect to a laboratory 
axis. The {'") is a Clebsch-Gordan coefficient and Cim{R) is the modified spherical har- 
monic function [SS]. Inserting Eq. ^ into Eq. ([T]) and after a number of transformations 
following Ref. [37] we find that the dispersion potential is the sum of the isotropic ?7disp,iso(-R) 
and anisotropic ?7disp,aniso(-R) potentials, which can be expressed in terms of the molecular 
dynamic polarizability tensor at imaginary frequency, a'^^^{iu), of molecule A and B. 
We find an isotropic f/disp,iso(-R) potential 



^ISO 



^iso (Jii 

f/disp,iso(i?) = -^6m^,m'a^Mb,m'b - ^('A,MA;tBMB\To^o{'2,2)\tA,M'^;iBM's) , (3) 
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where the rank-/ tensor operator Ti^rni{k,p) is defined by 

{tA,MA;iB,MB\Ti^mXk,p)\iA,M'^;iB,M's) = I Y.^lmi\kpqr) (4) 

X {JaMa I JAkM'^q) {JbMb \ JBpM'j,r) 

and the normahzation Ni{k,p) = {I0\kp00) (JAJAlJAkJA^) {JbJb\JbpJbO) , such that the 
operator (. . . \Tijni{k,p) \ • • •) is one for the stretched state \iA, Ma) = \iA, M'^) = \X, vJaJa) 
and Iib.Mb) = |^_B,M^) = \X,vJbJb)- Here (j3m,3|j2ji"^2"^i) is a Clebsch-Gordan coeffi- 
cient. Consequently, 

= - / du{iA, JA\a'^{iuj)\iA, JA){'i'B, JB\a^{i(^)\iB, Jb) 

TT Jo 

^6 22 = -jr~ / duj{iA,JA\^oi'^{ioj)\iA,JA){iB,JB\^a^{ioj)\iB,JB) 
' 457r Jo 

where, for each molecule A and B, a = {axx + oiyy + azz) /"i and Aa = azz — {oixx + oiyy) /2 are 
given in terms of diagonal x, and z components of the polarizability tensor at imaginary 
frequency. For atoms and diatoms the polarizability tensor is fully determined by a and 
Aa. 

The anisotropic van der Waals potential is defined by 

f/disp,aniso(i?) = ^E(-irC'2„.,(i?) {Cg^^" (^A, M^; , |T2,_^, (0, 2) K^, M^; z^, M^) 



nil 



+ QT(2A, Ma, IB, MBlT^^^rniC^, 0)\tA, M'^; tB, M'^) 

+C6,22°(m, Ma; tB, MB\T2,-mi (2, 2) |u, M'j,; tB, M'^)} 



and the van der Waals coefficients are 

1 /-"o 

TT JO 



Qol" = / du{iA, JA\a^iiUj)\iA, JA){iB, JbI^c,^ iiuJ)\^B, Jb) 

TT JO 

1 r'^ 

^620° = / dco{iA,JA\^a^ii(^)\iA,JA){iB,JB\a^iiuj)\iB,JB) 

TT JO 



TT Jo 

y^aniso '^/^iso 

^6,22 ~ Y^6,22 ■ 

In our framework the x, y, and z component of the diagonal dynamic polarizability at 
imaginary frequency is determined as 

{z,M\annM\t,M) = ^ Y: ,^ ^^^^2^i^,,^2 ^\^S\d.h\i,M)\' (5) 

where h is the unit vector along the n=x, y, and z direction, and i, M and S denote rovibra- 
tional wave functions of a single molecule, {S\d\i,M) are matrix elements of permanent or 
transition electronic dipole moments. Equation ([s]) includes a sum over the dipole transitions 
to the rovibrational levels within the ground-state potential as well as to the rovibrational 
levels of electronically-excited potentials. Contributions from scattering or continuum states 
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of the electronic potentials are also included. Finally, c is the speed of light and eo is the 
electric constant. 

In summary, we have found various contributions to the long-range van der Waals po- 
tential. They are characterized by their angular momentum dependence. In particular, the 
term proportional to Cg" is independent of magnetic projections and the relative orbital 
angular momentum between the molecules. The term proportional to Cq22 induces coupling 
between the magnetic projections of the two molecules without affecting the relative orbital 
angular momentum. Finally, the three Cg°^*^° contributions do cause mixing between mag- 
netic projections and relative orbital angular momentum. Note that the coefficients Cq22 
and Cg will be of the same order of magnitude. They both are proportional to a product 

of two Aa. Dispersion terms proportional to the spherical harmonic Cim{R) for / > 2 also 
exist. We do not consider them here as they are smaller. 

III. VAN DER WAALS COEFFICIENTS OF POLAR MOLECULES 

This section describes our results for the van der Waals coefficients and Cq^^^° between 
two polar KRb dimers and two RbCs dimers. These molecules are prepared in rovibrational 
states of the ground X^E+ potential and are of interest for on-going ultracold experiments [H 
El [23], SI] . We have used a non-relativistic ab initio version of the electronic structure multi- 
reference configuration interaction method ^39l HQ] to obtain potential energies, permanent 
and transition electric dipole moments of the KRb and RbCs molecules as a function of 
internuclear separation. In order to determine experimental observables we combine our 
electronic structure calculations with that of the calculation of rotational- vibrational energy 
levels to find vibrationally-averaged transition dipole moments. 

The isotropic Cq° coefficients between two molecules both in the J=0 rotational levels of 
the X ground state potential as a function of vibrational quantum number are shown in 
Fig. [Tj Comparison of the two panels show that for the deeply bound J=0 vibrational levels 
the isotropic coefficient for RbCs is almost an order of magnitude larger than for KRb. This 
is because for these levels the main contribution to the total value of Cq° is due to transitions 
within the ground state potential, which to good approximation is equal to df,/{6By) |36] . 
where dy and are the vibrationally averaged permanent dipole moment and rotational 
constant of level v, respectively. The permanent dipole moment of RbCs is twice as large as 
that of KRb, while is almost half the size [21 EHl 112] • The permanent dipole moment of 
both molecules rapidly decreases to zero with v and the X ^S"*" contribution to the van der 
Waals coefficient follows this trend. Consequently, for the highly excited vibrational levels 
the Cq° are solely determined by transitions to electronically excited states. In fact, for 
KRb the excited state contribution plays a dominant role for all v. 

Figure [T] also shows the J=l isotropic dispersion coefficients. These are independent of 
the projection quantum number M. In this case the contributions from the transitions from 
J=l to J=0 and 2 within the ground state potential have opposite sign and nearly cancel 
each other. As a result, the value of Cq° is close or equal to the value obtained by only 
including the excited state contributions. The dispersion coefficient is weakly dependent on 
V. The isotropic values for the v = vibrational level of the X^S"*" state of KRb and RbCs 
are given in Table |T| 

The J=0 and 1 dispersion coefficients are the same when they are solely determined by 
transitions to electronically excited states. This is because rotational energy splittings are 
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FIG. 1: Isotropic molecule- molecule van der Waals coefficients in atomic units for the J=0 and 1 
rotational levels of the X^S"*" ground state of RbCs (left panel) and KRb (right panel) as a function 
of vibrational quantum number. The curves labeled "J = (transitions within X)" correspond to 
an isotropic van der Waals coefficient where only contributions from transitions to ro-vibrational 
levels of the X^S"*" potential are included. 



negligible compared to electronic excitation energies and the rotational energy dependence 
of Es^ and Esg in Eq. [T] can be neglected. Alternatively, this independence follows from 
a simple atomistic model for the dispersion coefficient of weakly bound molecules. The 
molecular electronic wave functions are then to good approximation products of atomic 
electronic wave functions and the molecule-molecule dispersion coefficient reduces to a sum 
of atom-atom dispersion coefficients, which do not depend on the rotational state of the 
molecule. In fact, for RbCs+RbCs we find that C'i° 2C6(RbCs) + C6(Rb2) + C6(Cs2) = 
22868 a.u. for large v using the well-characterized Rb+Rb, Rb+Cs, and Cs+Cs dispersion 
coefficients |13]. This value is in good agreement with our results. A similar agreement is 
found for C^'^ of the interacting KRb molecules based on values of Cg from 03]. 

The collisions between two J=l molecules are anisotropic. Once molecules are prepared 
in the specific state characterized by M, the interaction energy can depend on the relative 
orientation of the two molecules. In fact, the anisotropic interaction can change the pro- 
jection quantum numbers. Figure [2] (left panel) shows our results for the anisotropic Cg "I" 
coefficient as a function of vibrational quantum number for J = 1 rotational sublevels, 
M = and M = ±1, of the X^S"*" ground state of KRb and RbCs. The anisotropy for 
RbCs changes sign as function of v due to competing contributions from the ground and 
excited states. In both cases Cg 02° goes to zero for highly excited vibrational levels. For 



these collision we have Cg 20° 



/^aniso 
'-'6,02 • 



Our predicted isotropic van der Waals Cq" coefficients for the interaction between KRb 
molecules have been used in theoretical models [23], IMI to successfully describe the loss 
rate constant /C observed in a JILA experiment as a function of temperature [23]. The 
anisotropic Cg"^^*^" coefficients can be used for describing collisions between molecules in 
non-zero J rotational states. Our estimate shows that an experimentally-accesible external 
electric field of 2 kV/cm will induce an ^ 50 MHz splitting between M = and M = ± 
1 components of J = 1 rotational state of KRb. The anisotropic interaction terms, for 
example, with coefficient Cggf" ~ 2500 a.u., will then cause a reorientation of the magnetic 
sublevels for separations less than 80 oq and will contribute to loss of molecules from the 
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FIG. 2: Anisotropic van der Waals coefficients Cg (left panel) and Cg "f" (right panel) in atomic 
units as a function of vibrational quantum number for two molecules in a J = 1 rotational levels of 
the X^S"*" ground state of KRb (solid lines) and RbCs (dashed lines). The coefficient is shown for 
both M = and 1. The coefficient for M = follows from the one for the stretched state M = 1 
by evaluation of the tensor Timi in Eq. |4[ 

trap. 

We also determined the coefficients Cq22 C'g "I °. For rotationless J = molecules both 
coefficient are zero. Right panel of Fig. ^ shows C'g 22° ^'^^ J = 1 as function of vibrational 
quantum number. These values are an order of magnitude smaller than the contributions 
discussed so far. The anisotropic coefficients for both molecules in the f = X^S"^ state are 
given in Table [Tj 

It was observed in Refs. [T9H23] that the lifetime of both Feshbach and deeply bound polar 
molecules decreases when they are in the presence of ultracold atoms. This suggests that 
collisions between atoms and molecules play a role in limiting the molecule lifetime. These 
losses are due to relaxation of rovibrational or hyperfine degrees of freedom that have enough 
kinetic energy to remove both molecule and atom from an external optical or magnetic trap. 
Another possible loss channel is an ultracold chemical reaction. To understand these losses, 
we determine the isotropic and anisotropic van der Waals coefficients for the interaction 
between a molecule in state |zmo1, Ma) and atom in state \iAt, Mb). Again following Ref. [37] , 
the isotropic coefficient is 

TT JO 

where a^°' and d^^ are the mean molecular and atomic polarizabilities, respectively. The 
anisotropic coefficient is 

^AtMol,amso = 1 /" rfo; (ij^^j | Aa^°' (iw) |iMol) («At | tt"^* (^w) lUt) , (7) 
TC Jo 

ml 1 r ^AtMolJso ^AtMoKaniso i ^AtMoLaniso A • r j. 

ihe values tor Cg 22 , ^6,02 ; ^J^cl Cg 22 are zero as Aa is zero tor an atom. 

Figure [3] shows the isotropic van der Waals coefficients for RbCs and KRb molecules with 
constituent atoms of Cs, Rb and Rb, K, respectively. It provides evidence that the van der 
Waals coefficients for both J=0 and 1 are nearly the same. This is due to the fact that 
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FIG. 3: Isotropic van der Waals coefficients in atomic units as a function of vibrational quantum 
number for the collision between a ground X^S^ state RbCs molecule and either a Rb or Cs 
atom (left panel) and KRb molecule and either Rb or K atom (right panel). Rotational J=0 (1) 
levels are shown as solid (dashed) lines, respectively. For KRb the J=0 and 1 Cq° coefficients are 
indistinguishable on the scale of a figure. 



the contribution from transitions within the ground state are negligible. The lack of these 
contributions effects the behavior of an anisotropic coefficient as well. As it is 

shown in Fig. 4, the values of the (j^^^"^'^"^^^'^ for the same projections M have the same sign. 
Table |l] lists the isotropic and anisotropic values for the v = vibrational level of the X^E"*" 
state of KRb and RbCs with atoms. 

Similar to the molecule-molecule dispersion potential, the van der Waals coefficients for 
one molecule and one atom are additive for a weakly-bound molecule. In other words, the Ce 
coefficient for the highly excited vibrational levels can be expressed as a sum of contributions 
from individual di-atoms. For example, (^^'^Cs+Rb.iso _^ C*g(RbCs) -|- C6(Rb2) = 10035 a.u. 
and Cgi^bcs+cs.iso _^ C'6(RbCs) + C6(Cs2) = 12514 a.u [13] for weakly bound RbCs molecules. 
The agreement between the asymptotic values and our numerical calculations is satisfactory. 
The same estimate can be done for KRb molecules interacting at long distances with K or 
Rb atoms. This lead to 0^^^-^^^'^'° C6(KRb) + ^(Rbs) = 8965 a.u. and Q^i^^+^.i^o ^ 
C6(KRb) + Cq{K2) = 8171 a.u., which agree within a few % with the numerical values of 
Fig.m 



IV. REACTIVE COLLISIONS OF ULTRACOLD POLAR MOLECULES 

The lifetime of ultracold molecules in an optical trap is determined by the inelastic colli- 
sional loss rates between these molecules. A collision can change the internal rovibrational 
state of the molecules as well induce reactions where the bonds between the atoms are 
rearranged. For example, for an ultracold KRb gas the reaction KRb+KRb — )■ K2+Rb2 
is allowed |23]- Once the molecules are prepared in their absolute ground state only the 
reactive process is present. 

This section describes our results using a quantum reflection model of the rate coefficient 
for rotationless J = '^°K^''Rb and ^^Rb^^'^Cs ground state molecules. It is based on a 
modification of the approach developed in Refs. [3ll |35]. We solve a single-channel radial 
Schrodinger equation for the inter-molecular separation with an isotropic van der Waals 
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FIG. 4: Anisotropic van der Waals coefficients Cg 20° atomic units as a function of vibrational 
quantum number for J = 1 rotational M = and M = ±1 sublevels of the X^S"*" ground state 
of the RbCs+Rb (solid lines) and RbCs+Cs (dashed lines) interacting systems (left panel) and 
KRb+Rb (solid lines) and KRb+K (dashed lines) interacting systems (right panel). 



potential for R > Rs^, where i?sr is a short-range separation to be defined below. Two 
short-range parameters describe the collisional wavefunction when the molecules are "close" 
together at R = Rsr- As we will show this boundary condition models the physics or reac- 
tivity for R < Rsr and "partially" inelastic events can occur. The short-range parameters 
can depend on the rovibrational state of the molecules. The isotropic van der Waals coef- 
ficients are taken from Table [Tj Recently, a similar model with a different treatment of the 
short-range interactions was used in Ref. [21] for the collisions between KRb molecules. 

We use the van der Waals potential to determine the validity of the model and es- 
timate characteristic length scales to distinguish between the dispersion and short-range 
interaction regions. The range of the isotropic van der Waals potential is defined as 
R < Rq = ^2/iCf°//i^, where /x is the reduced molecular mass [H]. For t> = 0, J = 
KRb and RbCs molecules this separation is ~ 300 and 500 oq, respectively. The short- 
range Rsr is defined by Cq"/R% ~ 2i?„, where is the rotational constant of the ground 
state molecule. For shorter separations the energy of the dispersion potential is much larger 



TABLE I: Van der Waals coefficients in atomic units for the interaction between two molecules 
in the = 0, J = and 1 rovibrational levels of the X ^S"*" potential as well as between such a 
molecule with an atom. The uncertainty in the coefficients is 5%. 
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{J,M) = (1,0) 


(J,M) = (1,±1) 


(J,M) = (1,0) 
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KRb + KRb 


16133 
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-11 
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142129 


16865 
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-157 


KRb + Rb 


7696 


7686 


1428 


-714 






KRb + K 


6905 


6896 


1278 


-639 






RbCs + Rb 


7326 


7274 


798 


-399 






RbCs + Cs 


8479 


8416 


929 


-465 
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than the rotational energy and coUisional interactions can mix rotational states. Such mix- 
ing falls outside the scope of our single channel description. We find values for Rg,. between 
50 oq and 80 Oq for both molecules. Hence i?sr ^ -Re- In fact, C^°/R^^ is much larger than 
collision energies of interest as well as the centrifugal potential h'^l{l + l)/(2/xi?gj,) between 
molecules. The quantum number / is the partial wave between the molecules or one molecule 
and one atom. 

Our model involves scattering of rotationless molecules in the potential —Cq° / + h'^l{l + 
l)/{2fiR'^) for R > Rsr- The contribution for partial wave / and projection m to the inelastic 
rate coefficient is given by 

Kl'^^{E)=v^^Y.\S^c.{E,lm)\' = v^^{l-\S.,{E,lm)\') , (8) 

where E = h'^k'^ / (2fi) is the collision energy and v is the relative velocity. The scattering 
S matrix elements Sia{E,lm) describe the transmission and reflection amplitude from the 
initial state i of colliding molecules and atoms to a flnal state a with either molecules in 
different rovibrational states or with a different bond. The sum over a excludes the initial 
state. Flux conservation or the unitarity of the S matrix allows us to rewrite the loss rate 
coefficient in terms of the diagonal S matrix element, Su. At ultracold temperatures, only a 
few partial waves / contribute as for higher I the centrifugal barrier prevents the molecules 
from approaching each other and Kj'^^ rapidly goes to zero. 

The scattering S matrix is calculated from the radial Schrodinger equation 



\ 2^idR^ R^ 2/i R^ 
at collision energy E with the boundary condition 

iJiUR) = A (e+^t^-'^/^l - (iUE)e'''^-^^^e-'^y--/'^) (9) 

at R = i?sr <^R6,y= {R/R(s)-^/2, and A is a normalization constant ^44^- The exponents 
g±«[y-7r/4] recognized as WKB solutions of a van der Waals potential at zero collision 

energy. The short range amplitude CimiE)e'^^^^"^^^^ and thus the real parameters Cim{E) and 
5ira{E) are not known a priori. They can only be determined from the chemical bonding 
between the three or four atoms. However, once they are known the boundary condition 
uniquely specifies the solution of the Schrodinger equation. Flux conservation requires that 
< Cim < 1 , where C«m = corresponds to the case where no flux is returned from the short 
range and = 1 correspond to the case where everything is reflected back. The phase Sim 
describes the relative phase shift of the flux that returns from R < i^gj.- In Ref. [2l] a model 
potential with an optical or imaginary term is used to describe the short range physics for 

R < i?sr- 

The van der Waals potential is the largest energy scale at R = Rsr and, consequently, 
we initially assume that both Cim{E) and Sim{E) are independent of collision energy E, 
partial wave /, and projection m. The validity of these approximations can only be tested 
by comparison with experimental data. We will do so with the data on KRb loss rate 
coefficients from Ref. [23] . 

In the limit R ^ oo the wavefunction approaches 

iJi{R) -> e-*^'^"^'"/') - Su{E, /m)e*('^^-''^/2) , (10) 
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FIG. 5: The total inelastic loss rate coefficient for a nuclear spin unpolarized (left panel) and 
polarized (right panel) sample of f = 0, J = ^"K^^Rb molecules in the X^S^ potential as a 
function of the short-range parameters C and 6 at a collision energy oi E/ks = 350 nK. The loss 
rate of the unpolarized case contains non-negligible contributions from s and p wave contributions. 
The polarized case only contains odd partial waves. 

where Sii{Elm) is the diagonal S matrix element required for Eq. ([s]). In order to determine 
the relationship between OmC^*'^''" and Sii{E,lm) we solve the Schrodinger equation ana- 
lytically in the limit Rsr — )■ using the solutions from [15]. Alternatively, the Schrodinger 
equation can be solved numerically. 

The left panel of Fig. |5] shows the total loss rate coefficient summed over / and m for 
collisions between fermionic ^"K^^Rb molecules in the v = 0, J = rovibrational state of 
the X^S"*" potential as a function of the short-range parameters ( and 5 at a collision energy 
of E/ks = 350 nK. For this figure we assume that C/m = C cind 6im = S are independent of / 
and m. This rate coefficient can be compared with experimental data on an "unpolarized" 
sample of this fermionic KRb molecule with at least two equally-populated nuclear-hyperfine 
states. The right panel shows the loss rate coefficient at E/ks = 350 nK for a "polarized" 
sample where all molecules are in the same nuclear hyperfine state and only odd partial wave 
scattering is allowed. In fact, for this collision energy only the p wave is significant. The p 
wave loss rate is much smaller than that for the s wave, since the centrifugal barrier is larger 
than the 350 nK collision energy and the p-wave loss rate is reduced. We use Cq° = 16133 
a.u. for both panels. 

The most striking observation of Fig. [5] is the appearance of resonances in the loss rate. 
The origin of this behavior is in the interference between the in- and out-going flux for 
R > Rsr- The maximum at {( ^ 0.8, 6 ~ O.Gtt) is due to s-wave collisions and that at 
{( ~ 0.95, 6 ~ 0.97r) due to p-wave collisions. In fact, the loss rate is signiflcantly larger 
than the rate at C = when no flux returns from R < Rs^- (For C = the loss rate is 
independent of 6.) The loss rate can also be much smaller than at ( = 0. In the limit C ~^ 1 
we have K^^^^ — )■ 0. For clarity, this limiting behavior is not shown in the surface plot. If we 
learn how to change the short range parameter we can envision that a significantly reduction 
of the loss rate is possible. 

Figure [6] shows a comparison of thermally-averaged rate coefficients for ground state 
^°K^''Rb molecules with the experimental observations of Ref. |23]. The thermal average 
assumes a Maxwell-Boltzmann distribution. In order to obtain the comparison we not only 
assumed that the dm and Sim are independent of partial wave but also independent of 
collision energy over several kT, where k is the Boltzmann constant. 
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FIG. 6: Bounds on the allowed values of the short-range parameters C and 6 from the experimental 
loss rates obtained in Ref. [23] for either a nuclear spin unpolarized (black curves) or polarized 
(red curves) sample of f = 0, J = '^'^K^^Rb molecules in the X^S"*" potential at a temperature 
of T = 250 nK. The full lines correspond to pairs of (C, S) where the experimental loss rate agrees 
with our thermally averaged loss rate. The dashed lines follow from the experimental uncertainty 
limits. In fact, at a temperature of 250 nK the measured unpolarized loss rate is /? = 1-9(4) • 10^^*^ 
cm'^/s, while that for the polarized case is /3 = 3.3(7) • 10~^^ cm^/s. 



The experimental rates are /3u = 1.9(4) ■ 10"^*^ cm^/s for the unpolarized case and /3p = 
3.3(7) ■ 10~^^ cm^/ s for the polarized case at a temperature of T = 250 nK. In the experiment 
/3 was extracted from the time-dependence of the molecule number density of each nuclear 
spin state. Consequently, our theoretical loss rate has been multiplied by a factor two 
for a polarized sample to account for the loss of two molecules per inelastic collision, and 
multiplied by one(l) for the unpolarized sample as one atom of each nuclear hyperfine state 
is lost. Ref. [23j verified that the rates for the polarized and unpolarized sample satisify 
their f3^ oc const, and (3p oc kT Wigner threshold limits as a function of temperature upto 
T ^ 1 /iK. 

The figure shows regions in the plane (C;^); where the theoretical rate coefficient agrees 
with the experimental rate coefficients within the uncertainties. One region is for the p- 
wave dominated polarized case and one for the s-wave dominated unpolarized case. The 
two regions do no overlap. Hence, our initial assumption that the short-range parameters 
are independent of partial wave I is not valid. The validity of the Wigner threshold limit over 
collision energies upto a few /iK suggests that the short-range parameters are nevertheless 
independent of energy over this range. In other words the behavior at i? < Rsr, where 
a multi-channel description is needed to describe the reactive rearrangement of the atoms, 
must be captured by more than two parameters. The uncertainty in the dispersion coefficient 
does not modify this conclusion. 

Figure [t] shows the loss rate coefficient for the collision between two bosonic ^''RbCs 
molecules in the vibrationally excited v = 1, J = level of the X^S"*" potential at a collision 
energy of E/k = 250 /iK and Cq° = 142540 a.u. The nuclear hyperfine state are equally 
populated and partial waves up to / = 8 are included. We again assume that the short- 
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FIG. 7: The total inelastic loss rate coefficient for an unpolarized sample of u = 1, J = RbCs 
molecules in the X^S+ potential as a function of the short-range parameters C and 5 at a collision 
energy of E/ks = 250 jiK. At this collision energy many partial "waves contribute. 

range parameters are independent of partial wave. In the figure maxima or resonances for 
/ = 5, 6, and 7 are observed. As can be found from a comparison of the binding energies 
of the f = vibrational levels the three RbCs, Rb2, and Cs2 dimer molecules, chemical 
reaction at ultracold temperature RbCs+RbCs— )■ Cs2+Rb2 is endothermic by ~ 40 cm~^ 
and not energetically allowed. For collisions between f = 1 RbCs molecules the reaction is 
exothermic. At a 250 /xK collision energy rate coefficients are as large as 2 ■ 10^^" cm^/s can 
be seen in Fig. [7| With typical densities in ultracold experiments around 10^^ 1/cm'^ to 10^^ 
1/cm^ lifetimes as short as a 10 ms could be observed. 



We have performed a theoretical study critical for the physical realization of highly con- 
trolled ultracold molecular systems. In particular, we focused on hetero-nuclear KRb and 
RbCs molecules used in ongoing ultracold experiments. Our detailed calculations of Cg 
dispersion coefficients for interactions between molecules or molecules with atoms helped 
the experimental [23] and theoretical [M| efforts to quantitatively describe quantum-state 
controlled chemical reactions. 

Analyses of the isotropic and anisotropic interaction potentials of ultracold polar 
molecules unveiled a significant contribution from dipole coupling to electronically excited 
states. This leads to a dramatic change in the interaction potentials as compared to previous 
estimates based solely on the permanent dipole moment. 

Our scattering calculations predict constructive and destructive interferences in the 
molecular scattering loss rates as a function of short-range parameters. Comparison to re- 
cent experimental measurements p3] shows that the initial assumption that the short-range 
parameters are independent of partial wave / is not valid. 



V. CONCLUSION 
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